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SENSITIVITY OF SHORT -PERIOD TRACKING DATA 


FROM A LUNAR SATELLITE TO THE 

LUNAR GRAVITATIONAL FIELD HARMONICS 

By William R. Wells 
Langley Research. Center 


SUMMARY 


A study has been made to determine the sensitivity of short-period tracking 
data from a lunar satellite to the zonal harmonics up to degree four and to the 
first two sectorial harmonics (which are even functions of the longitude) of 
the lunar gravitational potential. The sensitivity of the tracking data is 
indicated by the differences which result whenever the range and range-rate 
values (relative to the center of the earth) are computed with and without var- 
ious gravitational components present in the gravitational potential function. 

A parametric study of the effect of inclination and nodal position on these 
sensitivities is also presented. 


INTRODUCTION 


The use of artificial satellites to determine the external gravitational 
field of the moon is currently under investigation. The success of this tech- 
nique will depend, to a great extent, upon the sensitivity of the tracking data, 
that is, range and range-rate measurements of the orbiting lunar satellite to 
the lunar gravitational field, and upon the ease with which this field can be 
separated into its various components. Once this separation has been accom- 
plished, the effect of each component on the range and range rate can be 
accounted for individually by means of a harmonic analysis. An indication of 
the sensitivity of the tracking data to each of these components is indicated 
in the difference in the calculated values of range and range rate with and 
without a particular harmonic present in the calculation. 

The purpose of this paper is to present an analytical determination of the 
sensitivity of the range and range rate of lunar satellites to various compo- 
nents of the lunar gravitational field during several satellite orbital periods. 
A comparison of these sensitivities to the tracking noise level should provide 
a preliminary indication of the response of the tracking data measurements to 
the lunar gravitational field. 



SYMBOLS 


Keplerian elements 

Aa,Ae,Ai,AoD,Afi,£M perturbations in Keplerian elements 
bp, . . . ,bej coefficients defined by equation (B22) 


C J " 


^nm'^nm 
dq > • • 
D 
E 

y • • 

F 


g 


1 ' * ’ 


h 

hqj • • 


cos^v dv ( j = 1, 2, . . . ) 

coefficients of lunar-gravitational-potential harmonics 
. ,d^ coefficients defined by equation (B23) 

mean distance between centers of earth and moon, 384,402 km 
e c cent ri c anomaly , rad 

. ,f^ coefficients defined by equation (Bl8) 

row vector formed by partial derivatives of range with respect 
to Keplerian elements 

. ,g coefficients defined by equation (B15) 

row vector formed by partial derivatives of range rate with 
respect to Keplerian elements 

angular-momentum vector of satellite, km^/sec 

. ,h^ coefficients defined by equation (B2l) 


IlO'Iu'WlK) integrals defined by equations (B7c), (B10c), (Bl4c), 

and (B20c), respectively 


^l> m l*ki direction cosines of the position vector relative to 

x f ,y f , z '-axis system 

— > ■ — ) 

m 2 ^k2 direction cosines of the vector h X r relative to 

x T , y 1 , z ! -axi s system 


^ 3 ^ m 3 y k 3 


direction cosines of the angular -momentum vector relative to 
the x^y^z’-axis system 


n 


mean motion of lunar satellite, rad/sec 
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mean motion of moon about earth, 0.27 X 10 ”5 rad/sec 
associated Legendre function 

distance from center of moon to satellite, km 

position vector from center of moon to satellite, km 

disturbing function due to a nonspherical nonhomogeneous moon, 
km^/sec^ 

mean radius of moon, 1738.1 km 

nm component of R as defined by equation (B3) 
time, sec 

lunar gravitational potential function, km 2 / sec 2 
true anomaly, rad 

Cartesian coordinates with respect to an inertial coordinate 
system 

Cartesian coordinates with respect to a moon-fixed coordinate 
system 

coordinates of center of mass of moon with respect to 
x 1 ,y T ,z '-axis system 

column vector formed by perturbation in Keplerian elements 

angle between line joining earth and moon centers and line 
joining moon and satellite centers, rad 

vernal equinox 

longitude of satellite measured in equatorial plane of moon from 
mean earth-moon line, positive eastward, deg 

product of gravitational constant and mass of moon, 

4902.8 km^/sec 2 

range of satellite measured from center of earth, km 
perturbation in range, km 

variation in range due to variation in coefficient of gravita- 
tional harmonic, km 



II 


p range rate of satellite measured relative to center of earth, 

km/sec 

Ap perturbation in range rate, km/sec 

5Ap variation in range rate due to variation in coefficient of 

gravitational harmonic, km/sec 

$ latitude of satellite measured from lunar equator, positive 

northward, rad 

ft* longitude of ascending node measured from mean earth-moon line, 

positive eastward, rad 


Subscripts: 

m,n mth order and nth degree of harmonic 

o nominal value 

A dot over a symbol denotes differentiation with respect to time. 


GENERAL CONSIDERATIONS 


Analytical Formulation of Problem 

A lunar satellite will experience small perturbations in its orbital ele- 
ments due to the influence of the higher order harmonics of the lunar gravita- 
tional potential function. These disturbances will also cause variations in the 
tracking data measurements since the range and range-rate measurements can be 
related to the osculating or time-varying elements of the satellite orbit. The 
difference in range and range rate computed with and without a particular grav- 
itational harmonic will be defined as the sensitivity of the tracking data from 
a lunar satellite to that harmonic. 


For this analysis the expressions for the range and range rate of the lunar 
satellite, in terms of the osculating elements, will be given relative to the 
center of the earth. In addition, it will be assumed that the lunar equatorial 
plane and the earth-moon plane are coincident and that the moon is assumed to 
revolve about the earth in a circular orbit. It can be shown that the assump- 
tion of coincidence of lunar equatorial and earth-moon planes is well justified 
with respect to range and range-rate measurements if the ratio a/D is much 
less than unity, that is, for the case of close lunar satellites. 


The range of the satellite is given by (see fig. l) ; 

p(a,e,i,o>,fi,M) = ^D 2 + r 2 - 


( 1 ) 


4 



where 


Z]_ = cos g = cos (03 + v)cos ft* - cos i sin ft* sin(o) + v) 

An expression for the range rate of the satellite can he obtained from a 
direct differentiation, with respect to time, of the range as given by equa- 
tion (l), that is. 


p(a,e,i ,to,ft ,M) 


r(r - E^]_) - Dr ^l 
(d 2 + r 2 - 2rDZ 1 ) 1 / 2 


( 2 ) 


The perturbation in range and range rate due to a disturbance can be 
obtained analytically by expanding equations (l) and (2) about their undis- 
turbed values in a Taylor* s series 


P = P r 


+ 

Aa 

+ ^a 

Ae 

+ ^a 

Ai 

+ ^a 

5a 

o 

5e 

o 

d± 

o 

5u> 


Ar> 


+ 

dn 


An 
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P = Pr 



Aa + |£ 

Ae 

dp .. 

+ — — Aa 
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Si o 
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dp 
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(3b) 


It will be assumed in this analysis that, over a few orbital periods, the 
perturbations in the elements are sufficiently small so that the second-order 

terms, o[(Aa)^j, in equations (3) can be dropped. Thus, the expressions for 

the changes in range and range rate can be formulated in the following linear 
forms: 


1 

PR° 

11 
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a 

a 

II 

3 

(**a) 

Ap = p - p Q = G 0 Ax 
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where the vectors F, G, and Aa are defined as 


F = 


G = 



dp 

dp 

Bp 

dp 

dp 

de' 

di' 

do/ 

dft' 

dM 

dp 

dp 

dp 

dp 

dp' 

de' 

di' 

do/ 

dft' 

dM, 


(5a) 


(5b) 


Aa = 



(5c) 


The elements of the vectors F and G can be obtained analytically 
through use of equations (l) and (2). If e represents any of the elements, 
a, e, i, cc, ft, or M, then 


dp 

de 




(6a) 


dp _ 1 d(pp) 
de P de 



(6b) 


where 


d(pp) 

Be 





Expressions for the partial derivatives, 
in appendix: A. 





(6c) 


Sr Sl l 
de' de ' 




are given 


Equations ( 4 ) express the first-order perturbations in range and range 
rate for any general disturbance. These expressions will be used, in this 
analysis, to compute the perturbations in p and p due to the higher order 
harmonics of the lunar gravitational field. 
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The elements of Aa are obtained after a discussion of the lunar disturb- 
ance function. Each component of the moon T s gravitational field and its effect 
on the range and range rate of a lunar satellite can be accounted for individu- 
ally by an harmonic analysis. The lunar gravitational potential can be written 
in a form similar to that recommended for the earth potential (see ref. l) as 



1 + 


oo n 
n=l m=0 


cos m9 + S nm sin m9 )P nm (sin 0) 


(7) 


The associated Legendre function appearing in equation (7) is computed by 
the equation 


p nm( sin 90 


cos m / 2 0 V (-l) t (2n - 2t) J sln 11 " 111 " 2 ^ 
2 n t J ( n - m - 2t ) J ( n - t ) i 


( 8 ) 


where 



if n - m is even 


and 


k = 



if n - m is odd 


The potential function is defined in equation (7) so that the motion of 
the satellite relative to the moon is determined by the vector relation 


= VU 


(9) 


The first term of U causes elliptic motion of the satellite about the moon 
which can be described by a constant set of Keplerian elements. The remaining 
terms, which make up the function known as a disturbing function, cause higher 
order perturbations to the elliptic motion which can be formulated in terms of 
time-varying perturbations in the elements. These time-varying perturbations 
are determined from a solution of Lagrange T s planetary equations given in ref- 
erences 2 and 3 as 


da _ _2_ SR 
dt na cM 


(10a) 
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L. 



■niiiii 


de _ 1 - SR _ l/l - SR 
dt na^e SM na^e 


di _ cot i SR _ esc i SR 

dt " na2/TTli " na2 * Q 

dcp _ -cot 1 SR + ^ 1 - e 2 SR 
dt na 2 /l~- e 2 Bi na 2 e 


dfl _ esc i Sr 

dt 2 >/7“ 2 Si 

na^ V 1 - 

dM _ n (l - e 2 ) SR _ _2_ SR 
dt na 2 e Se na Sa 


(10b) 


(10c) 


(lOd) 


(lOe) 


(lOf) 


The disturbing function R is defined as 


R S U - d 
r 

co n 

- £ (jfrJ (C nm cos m9 + Snm sin m0^P nm (sin (ft ) (ll) 

n=l m=0 

The function R can be expressed in terms of the angular orbital elements and 
the true anomaly of the orbit, by use of the following relations determined 
from figure 1 


sin 0 = sin i sin(o> + v) 


cos (j) cos 0 = cos 2' cos(o) + v) - cos i sin fl 1 


cos <f> sin 9 = sin fi ' cos(cd + v) + cos i cos Q ' 


sin(o) + v) 
sin(cD + v) 


(12a) 

(12b) 

(12c) 
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The perturbations in p and p 
^e.g., the harmonic with coefficient 


due to an individual harmonic of R 
Cnm) from equation (4) are 


Ap(Cnm) = F 0 (15a) 

Ap(Cnm) = G 0 (l3h) 


Values of Aoc,( Cnm) for various values of n and m are given in appendix B. 
Variations of these perturbations due to variations in the coefficient C nTn 
are given as 


Sfie(Cnm) - F 0 8C nD (14a) 

«p(Cnm) - So Mam (l*«>) 

ou nm 


Equations (l4) follow directly from equations ( 13 ) since it has been assumed 
that the perturbations in the elements are linear in the coefficients of the 
gravitational harmonics. 

The variations in Ap and Ap given by equations (l4) have been previ- 
ously defined as the sensitivity of the tracking data to the lunar gravita- 
tional harmonics. The determination of their behavior with time is the purpose 
of this analysis. 


ANALYSIS 0E RANGE AND RANGE-RATE SENSITIVITIES 


In the following analysis, a time history of the range and range -rate sen- 
sitivities will be considered for variations in the coefficients C^q, Cpp, 

^ 22 .> Cjo* anc ^ c 4o* The magnitudes of these variations are taken to 

correspond to the smallest values that cause variations in Ap and Ap which 
lie outside the assumed noise level at some time during the tracking phase. 

The orbit used for this analysis has pericentron and apocentron altitudes of 
46 and 185 O kilometers, respectively. The value of the argument of pericentron 
is taken as 0° and is assumed to be constant for the short times considered in 
the analysis. However, a parametric study is performed on the effects of 
inclination and nodal positions on the sensitivity of the tracking data. 

The time dependence of the range-rate and range sensitivities is given in 
figures 2 and 3 for three consecutive orbits. The tracking noise level is 
indicated in these figures as 15 meters in range and 0.002 meter per second in 
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range rate. These are values of the accuracy in range and range-rate measure- 
ments currently "believed feasible. 

Illustrated in figures 2(a) , 2(b), 3(a) > and 3(^) is the periodic behavior 
of the range -rate and range variations due to variations in the coefficients 
Cio and Cpp of 10~5. Since these two coefficients can be related to the 
location of the center of mass of the moon relative to the origin of an assumed 
coordinate system (x* = f' = % s ll* 2 1 = R^io)' their variations rep- 

resent an uncertainty in the location of the moon T s center of mass. Values of 
c 10 an ^ On of 10"5 correspond to an uncertainty of about 17 meters in the 

location of the moon’s center of mass along the x 1 - and z f -axes. The results 
of figures 2(a) and 2(b) indicate that values of C-^q and of the order 

of 10~5 may be detectable in the range-rate data during the first three orbits 
of tracking. Figure 3(a),? however, indicates that the variation in range does 
not exceed the assumed noise level for this variation in C]_q. 

The sensitivity of tracking data measurements to variations in C 2 q is 
given in figures 2(c) and 3 ( 0 ). The variation in range rate is above the noise 
level during the first orbit for a variation in C 20 of 10"^. The range vari- 
ation becomes greater than the noise level during the second orbit for this 
variation in C 2 q. The secular effect (linear change with time) of the second 

zonal harmonic is evident in the second and third orbits. 

The sensitivity of the tracking data measurements to a variation in the 
coefficient C 22 of 10"^ is given in figures 2(d) and 3(d). The variations in 
both range and range rate are beyond the noise level in the second orbit and 
continue to grow secularly with subsequent orbits. These results indicate that 
the radar measurement should be highly sensitive to the effects of the second 
sectorial harmonic of the lunar gravitational field. 

The sensitivity of the tracking data measurements to a variation in C^q 

of 10"5 is given in figures 2(e) and 3(e). These results indicate the range 
and range-rate measurements should be fairly sensitive to the effects of the 
third zonal harmonic during the first orbit of tracking. This sensitivity 
increases with additional tracking time as evidenced by the long-period varia- 
tions (variations with angular frequency co) in the range and range-rate values 
in the second and third orbits. 

The sensitivity of tracking data measurements to a variation in C^q 

of 10 “6 is given in figures 2(f) and 3(f)- The variation in and range rate, 
due to this harmonic, exceed the noise level after one orbit of tracking; the 
variation In range exceeds the noise level only after the second orbit. The 
long-period and secular effects of the fourth zonal harmonic on these values 
are evident in the second and third orbits. 

The effect of varying nodal positions on the tracking data sensitivity 
during one orbital period is shown in figures k and 5* The peak magnitudes of 
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the range and range-rate variations change little, with the exception of those 
due to the second sectorial harmonic, with nodal positions during one orbit of 
tracking. Figures 4(d) and 5(h) indicate the sensitivity of the tracking data 
to the second sectorial harmonic to be highly dependent on the initial nodal 
position of the satellite orbit. 

In a determination of gravitational constants, the ability to change the 
inclination of the satellite orbit plane may aid in the separation of highly 
correlated coefficients. This separation could be accomplished by the use of 
more than one satellite. The effect of varying the inclination on the tracking 
data sensitivity is shown in figures 6 and 7 for one orbital period. With the 
exception of the odd zonal harmonics represented by C^q and C^q, the sensi- 
tivities are slowly varying functions of the inclination. 

The sensitivity of the tracking data to the odd zonal harmonics is approx- 
imately proportional to sin i as indicated in figures 6 (a), 6 (e), 7 ( a ) > 
and 7 (e). 


CONCLUDING REMARKS 


A study of the sensitivity of short-period tracking data from a lunar 
satellite to the harmonics of the lunar gravitational field with coefficients 
Cio, Cpi, C 20 j c 22 j and ClfO has been performed. It was shown that 

the range and range-rate measurements, relative to the center of the earth, are 
sufficiently sensitive to the first, third, and fourth zonal harmonics so that 

variations in their coefficients of the order 10“5 cause variations beyond the 
noise level during one orbit of tracking. In the case of the second-order 

zonal and sectorial harmonics, a variation of the order of 10 will cause a 
variation in range and range rate beyond the noise level. In all cases, the 
sensitivities of the tracking data are greatly amplified after three full 
orbits of tracking. This amplification resulted from the secular and long- 
period effects of the harmonics. 

It was shown that the effect of varying the nodal positions changed the 
maximum sensitivity of range and range-rate measurements to the gravitational 
harmonics only slightly, with the exception of the second sectorial harmonic. 
The effect of varying the inclination was to cause the sensitivities to the 
gravitational harmonics to vary slowly with the exception of the first- and 
third-order zonal harmonics, which vary approximately as the sine of the 
inclination. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 23, 1965 . 
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APPENDIX A 


EVALUATION OF THE PARTIAL DERIVATIVES OCCURRING IN F AND G 


cDr 

The partial derivatives — > 

oe 


needed for the evalua- 


^1 dr . 

ST* a? ” d ST 

tion of the vectors F and G defined by equations (5a) and (5b) can be 
obtained analytically from the expression for two-body elliptical motion 
written in terms of the osculating elements. The two-body results needed are 
taken from reference 2. They are 


(l - e cos E) 

(Ala) 

cos E - e 

(Alb) 

1 - e cos E 


V 1 - e^sin E 

(Ale) 

1 - e cos E 

E - e sin E 

(Aid) 


The following direction cosines of the vectors r , h X r, and h rela- 
tive to the x',y ! ,z T -axis system are given for convenience since they will 
occur frequently in the analysis: 


l i 

= cos(o) 

4* v)cos Q 1 - 

cos i 

sin ft ' 

sin(flD + 

v) 

(A2a) 

m l 

= cos(o) 

+ v)sin n 1 + 

cos i 

cos ft ' 

sin(tt> + 

v) 

(A2b) 



k^ = sin i 

sin(cn 

+ v) 



(A2c) 

z 2 = 

-cos ft ' 

sin(oo + v) - 

cos i 

. sin ft 

1 cos(<d + 

v) 

(A5a) 

m 2 

cos ft ' 

cos i cos(o) 4 

v) - 

sin ft ' 

sin (cd + 

v) 

(A3b) 



k 2 = sin i 

cos(o) 

+ v) 



(A3c) 
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l j = sin i sin Si 1 


m = -cos Si T sin i 


k^ = cos i 


Evaluation of 


^ and 
Se 


S Hi 
Se 


From equation (Al) 


Sr _ r 
Sa a 


— = - S_(cos E - e) = -a cos v 

Se 


Sr _ Sr_ Sr 

Si Scd Sfi 
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— e sin E 
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(AUa) 

(A4b) 

(A4c) 


(A5a) 

(A5b) 

(A5c) 

(A5d) 


From equations (Al) and (A2) 


where 


s l 1 


Sa 


= 0 


hl l _ , Sv 
Se 2 de 


Sv 

Se 


a sin E 

r \J 1 - e 2 


II + 



+ e cos v) 

1 - e 2 


( A6a) 


(A6b) 
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Bln 

= k, sin ' 

Si ^ 


(a6c) 
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-m 


1 


(A6d) 
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d*l _ ? Sv 
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Evaluation of 
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From equations (Al) and (A2) 


r = l&i sin E 
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(A7b) 
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SZ 
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1 _ 


= n-^k^ cos A ' + vkg sin fi ' 
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Sen 12 1 


Sz L _ S i 1 
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si]_ 

Sm 
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(A8i) 
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where 



II 


APPENDIX B 


EVALUATION OF Ax 


The elements of the vector 


cAx 

Sc 


can he obtained once a solution for Ax 


nm 


from the Lagrange planetary equations is accomplished. One form of the solu- 
tion to these equations is given in reference k. For the limited number of 
harmonics treated in this analysis ^i.e., the harmonics with coefficients C]_q, 

Cqq, C20 j ^22; CJOj and Ci K)) it was found to be convenient, as well as 
instructive, to develop expressions for the components of Ax as follows. The 
integrated form of Lagrange's equations for an individual component of the lunar 
potential will be written as: 


Aa 


"2 
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R nm 
- tr 


(Bla) 


Ae = 


Rr 
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o p nm 
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cot i Ri P 2 

— - . - + — R 
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2 e 
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-it. 


(Bid) 
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— t 


R-i CSC i 


[na ^ yj 1 - e^j 


(Ble) 
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A _ (l - e2) 
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( Blf ) 
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where 


and 



(B2a) 


(B2b) 


R i = 



d R nm 

Si 


dt 


(B2c) 


% = 



SHnm 

Sod 


dt 


(B2d) 


% = 



d R nm 


dt 


(B2e) 



cos m9 P nTT1 ( sin $) 


(B3) 


Note that is the n^m^ R component of the even part of the disturbance 

function. 

The use of equations (Bl) and (B2) will be illustrated for the case of the 
perturbations caused by Rpo* 

Perturbations Due to R^q 

The disturbance function R-^q is given by equation (B3) as 


R 10 


p r m 

n 2 




sin i sin(u> + v) 


(B4) 
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The partial derivatives to be substituted into equation (B2) are computed 
as follows: 


SR 


10 


Sa 



sin i sin(cn + v) 


(B5a) 


d R 10 

Se 



— sin i cos(cd + v) 
Se 


M-% „ ( a\ 2 . . „/a\ . / , % , (2 + e cos v)sin v cos(o> + v) 

= — — C-inl-l sin 1 2[— jcos v sm(a> + v) + i i 
a 2 U\r) |_W p _ e 2 


(B5b) 


Sr io 

Si 


^10 

Sen 



cos i sin(o> + v) 


sin i cos (od + v) 


(B5c) 


(B5d) 


^ R 10 _ d R ip 

SD Sf2 ' 


0 


(B5e) 


These results can now be integrated with respect to time to obtain the 
proper values of R a , R e , R^, R^, and Rq as expressed by equation (B2). 

Since v is the independent variable in these partial derivatives, it is 
advantageous to integrate with respect to v rather than t. This integration 
can be accomplished by using the two-body result 

I - n/TT^) 2 (*> 
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Then 


■p = 1 r (rf 5R 10 

a "n /ir^2 J U) Sa 


10 dv = -. nJ ^ 10 ■ sin i cos(a) + v) (B7a) 


\l 1 - e2 


R= 


_ 1 /"* ( r ^10 

n\/ 1 - e2 J W 


y/l - e 2 ^ V a / 
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i (1 + — cos v] 


“ (l-e 2 ) 5/2 

+ ^2v + eCp - 2C 2 - eC-^sin <nj 


cos^v 


cos o> 


(B7b) 


where 


L 10 = 


2na%C 10 
(1 - 


sin i 


^C 2 + eC-^sin to - + 2. cos vj 


cos^v cos CO 


(B7c) 


R-! = - mCi ° ■ cos i cos(cd + v) 


1 " /TT^2 


(B7d) 


naRi^rC-j q 

= 7 sin i sin(co + y) 


\J 1 - e 2 


(B7e) 


% = 0 


(B7f) 


The values of (j = 1, 2, . . .) appearing in equations (B7b) and (B7c) 


are defined as 


Cj = J' cos^v dv (j = 1, 2, . . .) (B8a) 
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These expressions can he computed, once and Cq are determined, by the 

recursive relation 


Ca~\ tan v co 


sJv + — 


C j-2 


(j “ 3> ^ } • • •) (B8b) 


Substitution of the results of equations (B7) into equations (Bl) provide 
the desired expressions for the elements of Ax. 


Perturbations Due to R- 


11 


The disturbance function R-j_-j_ is 


R 11 = ^ [f OS cos ( CD + v ) “ cos i sin H T sin(o3 + v)J 


Proceeding as before, the following results are obtained: 


(B9) 


*a = " 


2nR M c n 
\/ 1 - e 2 


jcos n' sin(cJD + v) + cos i sin fi’ cos (cd + v)J (BlOa) 


R e = I 11 " 


sa%Cu 


(1 - 


- (cos n 1 sin cd + cos i sin fi ’ ) ^1 + ^ cos 


cos^v 


+ (cos cos cd - cos i sin ft’ sin cd)^2v + eC^ - 2C 2 - eC^ 


(BlOb) 


where 


L 11 


2naR M C 11 


|(cos Q' cos cjd - cos i sin O' sin cd) ^C 2 + eC^ 


(l - e^ 2 

+ (cos 1 sin cd + cos i sin O' cos CJD )^~ + j cos v^cos^vj (BlOcJ 
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naRurCi 1 . . 

R,- = - ■ — sin i sin £2 * cos(oo + v) 

\fl - e 2 


(BlOd) 


R^ = - n 3* . U _ cos £2' cos(aa + v) + cos i sin £2' sin(co + v)J (BlOe) 

\j 1 - e 2 


r^ = _ |sjn £2' sin(a) + v) - cos i sin £2' cos(od + v)J (BlOf) 


yl - e 2 


Perturbations Due to R 20 

The perturbations due to R ^ are taken from existing expressions as 

given in reference 3* The only difference is that the present analysis includes 
the secular term in the disturbing function which is not included in the analy- 
sis of reference 3* The disturbing function for the present analysis is 


R 20 “ " 


^%^20 /a r 2 

4a3 Vi L 


- 3 


sin^i 


+ 3 


• 2- 

sxn^i 


cos 


2(to + v)] 


(Bll) 


while that used in reference 3 is 


R 20 


_1_ 

2jt 



dM 


or 


R 20 


+ ^C20 R M 2 / g 

4a3(i - e 2 f /2 


3 sin 2 i 


) 


For this case, the expressions given in equations (B2) are not considered 
since the perturbations can be written directly from reference 3* The pertur- 
bations in the elements due to the second zonal harmonic are 
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Aa = - 


3C 2 qRm J2 i 

2a 

3 


^ - 1 sln2i ) (I) 3 - 1 1 - e2 ) 


■3/2 


+ J* sin^i cos 2 (cd + v) 


For purposes of illustration is computed from equation (B12a) as 

SC 20 


cJAa 


c)C 


20 




v-3/2 


Ae = 


+ sin^i cos 2 (cjd + v) 

. 3Cgo% g (i ^ . 1 sln2l ) (ij . (! 

3 V 

+ ^r) cos 2(tu + v) > 

J v o 

+ V C 20 sin2 iL s 2(aJ + v ) + e cos(v + 2 co) 

4a2 e (l - e 2)L 


)•*] 


\ e cos(5v + 2 cd) 

5 _ 


Ai = - 


3 c 20 R M 2sin 21 


8a' 


% - * 2 f 


c< 


cos 2 (cd + v) + e cos(2cd + v) + — cos(an + 3v) 


) 


(B12a) 


(B12b) 


(B12c) 

v 

v o 

(B12d) 
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zNM = 


tv - - ?Cg0 ^ 2 ■ ■ - f( 1 - 2. sin 2 i)( 

4a 2 (l - e 2 ) 2 11 1 


v + e sin v) 


(l - - sin 2 i)["- 

-(l - i- e 2 ^sin 

v 2 'L e 

V 4 ) 

ifi- sin 2 i + ( — 

- ^5. sin 2 !^ 2 ] 

e[4 \2 

16 J J 

— sin 2 i sin(v 

- as) _ i/i - 

16 

2 \ 

i|X sin 2 i - H 

^1 - sin 2 iR e 

VO 

CV1 
1 — l 

<u 

^ 8 J 


— sin 2v + — sin 3v 

2 12 J 


J 


sin(3v + as) 


V 


+ — sin 2 i sin(4v + as) + — sin 2 i sin(5v + 2s) > (B12e) 

8 16 | 

J v o 


Aft = 


3C 20 %‘ 


2a‘ 


% c ° s 1 r 

^ - - 2 f L 


+ e sin v - i- sin 2(cs + v) 


- sin(v + as) - 2. sin(3v + as) 

2 6 _J 


(B12f) 


2a‘ 


D fl . I sin 2 i\i 

2(1 - e ^ 3 2 > 


5 ^ 20 % 


R 


2\ 2 
1 - Jsin v + % sin 2v + — sin 3v 

4 ) 2 12 


t G 2a 2 e(l - e 2 ) 
+ sin 2 i 


,3/ 2 


K 1 - 


sin 2 i^ x 


i-^1 + ~ e^sin(v + 2s) 


— sin(v - 2co) - — [l 

l6 


e2 


— sin(5v + 2s) 
Id 


0} 


1_ 

12 ' 

v 


|g-Jsin(3v + 2s) 


- e sin(4v + 2s) 

8 


(B12g) 
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22 


The disturbance function R 22 is 


R 22 = 


^22 /a r / \ 

2^1 gT^r) |fl “ S 2 + + cos v ) sin v cos v 

+ - g 2 )cos v + 2gg cos 2 v + 2egg cos-^vj 

The remaining expressions to be substituted into equations (Bl) are 


(B13) 


R a = - 


9 R m nC 22 r 2 


p e i)3/2t 3 683 COS3V " §3 + ( §1 " g2 ) V 

+ e ( s l " g 2) C l + 2g 2 C 2 + 2 e S 2 C 3 ] 




3n%*C 22 , _ 

e " 22 ' 2 g 5 


h gj cos v |a + ^ e cos V + j^e 2 - ^Jc° s2v ~ 
- 2 e 2 cosVj + ag 2 Jac, L + 3eCg + (e 2 - a)Cj - 3eCi, - e 2 C 5 j 


(1 - c 2 )= 


(Bl4a) 


— e cos^v 
2 


( Bl4b ) 


where 


X 22 = J 3^R 22 cos v dt 

9nR M 2 C 22 loo\ 

T/T <" 2g-,( “ + - e cos v + — e^cos 2 v]cos5v 
2(1 - e 2 ) 5/2 l A3 2 5 ) 

+ (gq - S 2 )Cq + 2e( gl - g 2 jC 2 + £e 2 ( g;L - g 2 ) + 2g 2 jc, 

+ 4eg 2 C^ + 2e 2 g 2 C 5 J 


(Bl4c) 
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Ri = . 3 ff“ ° 22 |7sin 2i cos 22' - gjAv 

2(l-e2) 5/2 L V ^ 

+ e^sin 2i cos 22’ - g^C-^ + 2g ij C 2 + 2eg if C^ - 


p 2 ^ 

5 ^ cos^v - — eg^ cos-'v 


(Bl4d) 


*Vd = _ g 2 ^l + | e cos v)cos 2 v + g^v + eC L - 2C 2 - 2eC 3 ) 


3n%2c 22 T 


Rn = - — - ■ ■ ^ , I- g v (l + % e cos v)cos 2 v 

(l-e2) 5/2 L n 5 


,|l + ^ e cos 

+ ^sin 2 i sin 22' - gg^v + e Cp^ + 2gg^C 2 + eC^jj 

In expressions (B13) to (Bl4f) the following definitions were employed 

= sin^i cos 22 ' 
g 2 = cos 2xi cos 22' (l + cos 2 i) -2s 
g 3 = -sin 2x cos 22 '(l + cos 2 i) - 2 

g^ = 2 sin i sin 22 ' sin 2x - sin 2i cos 22 ' cos 2x 

g 5 


(Bl4e) 


(Bl4f) 


g 7 



(B15a) 

22' cos i 

(B15b) 

. 22* cos i 

(B15c) 

cos 2x 

(B15d) 

cos 2x 

(B15e) 

2' sin 2x 

(B15f) 

2 ' sin 2x 

(B15g) 
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Perturbations Due to R 


30 


The disturbance function R^q is 


R 


30 


= ^ R ^50 ^ ^5 sin 2 i - ^sin(o) + v) - 5 sin 2 i sin 3(oo + v)j sin i 


(B1 6) 


The remaining expressions to be substituted into equations (Bl) are 


R„ 


= - — ^ — - - — [^e^sin^i cos 3^ cos^v + lOe sin^i cos 3^ cos^v 


2a' 


\l - e2) 


5/2 L 


+ ~^20 sin^i cos 3^ - e^fg^cos^y - efg cos^v - fg 
+ 5f]_C^ + 6ef^Cg + “ 20 sin 2 i sin 3^0^ 

- 40eCij. sin2i sin 3^ - 20e2c^ sin^i sin 


cos v 


(B17a) 


R= = 


— — — — — f^Oe^sin^i cos 3^° cos^v + 80e2sin2i cos 3 <*> cos^v 

8a(l - e2) 7 / 2 l 

+ £l08e sin2i cos 3^ - ^ e^4f 2 + 3f^jjcos5y + j^O sin^i cos 3 00 

- 3e2^fg + f ^ jj cos^v - e^4fg + 5f^^cos5y - ^2fg + 3f’^cos2y - 6f^v 

- 15ef 1 C 1 + 6c 2 (t k + 2f ± - 2e2f L ^ + 3e ^>f k + f x (l2 - e 2 )^ 

+ 2^ 100 sin 2 i sin 3^> + ^e 2 ^ + 3f ] _jjci J . + 3eje 2 ^f^ + kf^ 

- 180 sin2i sin 3“jo^ - 480e2sin^i sin 3^Cg - l40e3Cy sin 2 i sin 3^ 


(B17b) 
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R-? = 


- — R a cot 
4 a 


nRvr C^n sin 2 i cos i 0 „ ^ ^ , i, 

i + — ii — ^ = — { 4e^cos 3^ cos^v + lOe cos 3^° cos t 

fca(l - e 2f/ 2 I 

+ ^jh cos 3 o) - e^(cos 30 + 3 cos 0))j cos 5 v - 5 e(cos 30) + 3 cos o)cos^v 

- 5(cos 3^ + 3 cos cs)cos v + 15 (sin 3°° + sin CJD )C]_ + 30 e(sin 3°^ + sin u>)C 2 

+ 5j3e 2 (sin 3co + sin cjd) - 4 sin 3<nJcj - 40eCj^ sin 3u> - 20e 2 C^ sin 3°^ 


(B17c) 


R, n = n %l ^ C j 30 sin r i 2 e 2 sin 2 i sin 300 cos^v - 30 e sin 2 i sin 3 n> cosV 
8 a(l-e 2 ) 5 / 2L 

+ ^f^e 2 - 20 sin 2 i sin Jcn'jcos^v + 3 f]_ e cos 2 v + 3 f]_ cos v + 3 f^C 
+ 6 ef^C 2 + 3 ^ e ^fj - 20 sin 2 i cos 3 a ^Cj 
- 12060 ^ sin 2 i cos 3 ^ - 60 e 2 C^ sin 2 i cos 3 °^ 


3 1 


(B17d) 


% = 0 


(B17e) 


In equations (Bl6) to (B17e) the following definitions were employed: 


fp = ^5 sin 2 i - 4^ sin cjd + 5 sin 2 i sin 3^ (Bl8a) 
fg = 3(5 sin 2 i - 4^cos 0) + 5 sin 2 i cos 3^ (Bl8b) 
fj = (5 sin 2 i - 4^ cos cn + 15 sin 2 i cos 3“ (Bl8c) 
f ^ = (5 sin 2 i - 4) sin cu + 25 sin 2 i sin 3o> (Bl8d) 
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Perturbations Due to R^q 
T he disturbance function R^ is: 


R 4o = 


55 m-% c 4 o/a \5 


5 + ^2 cos ^ v + r 3 cos^v + h 4 cos^v sin v + sin v cos vj 

(B19) 


The remaining expressions to be substituted into equations (Bl) are: 


R q = - 


175nR M ^C 4o 


8a5(l - e 2) 7/2 


i e^h), cos^v - i e^hj, cos^v 


lcos + v - ehcr cos?v - - h R cos^v 

P 2 P 


- i e^e 2 h^ + Jh^cos^v - j^h 4 + 3e 2 h^< 

+ h^v + 3eh^Cn_ + ^hg + 3e 2 h-^C2 + e^hg + e 2 h^Cj + + 3e 2 hg^C 4 

+ e^3hj + e 2 h 2 ^C 5 + 3e 2 h 5 Cg + e^hjCyj (B20a) 


R~ = 


I 4o 


33nR M 4 C 4o 

- 2 'l - e 2 f/ 2 




b^e^cos^v - b-re^cos&v 


-* b^e^cos^v - i e^(e 2 b 0 + 9b^cos^v 

8 a 2 (l - af'‘ ^ ' ' 8 7 \ 2 

- | e^5e2b 2 + 7b^cos 6 v + “ 9e2b 2 - 2b-^cos5v 

- ^ e^e^b^ + Tb^cosV - j^ 9 e ^]_ + 2 b 2 ^cos 3 v - i eb-^ cos^v - 2 b^ cos v 
+ 2b 4 C 1 + 7eb 4 C 2 + (9e2b 4 + 2b 5 - 2b 4 ^C 3 + e [j(b 5 - b^ + 5e 2 b^j C 4 

+ |b 4 e^ + 9e 2 (b 5 - b 4 ) - 2b^jc 5 + e |5e 2 (b 5 - b 4 ) - Jb^C 6 
+ e 2 |9b 5 + e 2 (b 5 - b^jjc ? - ^e^Cg - b^C^j (B20b) 
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where 


x ko _ 


175nR^ Ci^q 


— e^K cos^v - — e^h k cos®v - e 2 (e 2 h,- + 6hh) 

Sa^l-e^L 9 24 7 1 5 V 

_ £ e^e 2 h^ + h^cos^v - ^6e 2 h^ + h^cos^v - eh^ cosV - j h^ 

+ h-jC-L + beh 1 C 2 + (h 2 + 6e 2 h-^jC^ + ^eC^hg + e 2 h 1 'j 
+ ^h.3 "*" 6e 2 h 2 "*" "*" 

+ e 2 ^ 6 h^ + e 2 h 2 ^Cy + 4 e 5 hjCg + e^h^C^J 


cos 'v 


cos-^v 


(B20c) 


R i 


35 nR M i 4 ’Ci J .Q 


- i dce^cos^v - — dc-e 2 cos^v - — efe 2 dh + 3d£\cos5v 

7 5 2 5 5 V 4 5 / 


8a 2 (l - e2) 7 / 2 L 
- i^d^ + 3e 2 d^cos^v - ed^ cos^v - i- d^ cos 2 v + d^v + 3ed-j_C-]_ 

+ ^d 2 + 3 e 2 d 1 ^C 2 + e^ 3 d 2 + e 2 d^C 5 + + 3 e 2 d 2 ^ 

+ e^3<3-j + e 2 d 2 ^C^ + 3 e2 djCg + e5djCyJ (B20d) 

R. = _ — 55nlfy[ Ckp [ i 3, cos 7 v - — h C re 2 cos^y - — e^e^bj, + 3bcAcos5v 

8a 2 (x-e 2 f/ 2 L 7 5 2 5 5 l * 5j 

_ + 3e 2 b^cosV - eb^ cos5y - £ b^ cos 2 v + b-jV + 3eb 1 C 1 

+ ^b 2 + 3e 2 b^ C 2 + e^3b 2 + + ^>3 + 5e 2 b 2 ^C^ 

+ e^ 3 b 5 + + 5e2b 3 C 6 + e 5 b 3 C 7 j ( B20e ) 


% “ 0 


(B20f) 
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In equations (B19) to (B20f) the following definitions were used: 


2- sin 2 i + 2: sin^i 

7 8 


cos 2a> sin2i[2. 

\7 


— sin 2 il + — sin^i cos to 
2 / 8 


(B21a) 


hg = 2(1 - sin 2 i^ cos an sin 2 i - cos to sin^i (B21h) 


= cos to sin^i 


(B21c) 


hi, = -sin to sin^i 


(B21d) 


he = - 2 [ — - — sin 2 i) sin an sin 2 i + — 
5 \7 2 ) 2 


— sin 2 i] sin an sin 2 i + — sin to sin^i (B21e) 


b-, = -2 (— - — sin 2 i\ sin 2cn sin 2 i + — sin^i sin to 


(B22a) 


4 sin 2 i ( 2 . - i. sin 2 i| sin 2cn - sin to sin 2 i 
\7 2 


(B22h) 


bj = 4 sin^i sin to 


(B22c) 


h^ = 2 2(1 - sin 2 i^cos an - sin 2 i cos to sin 2 i 


(B22d) 


he = 4 sin^i cos to 


(B22e) 


sin i cos i - — + 2. sin 2 i - 2(2. - sin 2 i | cos 2u> + — sin 2 i cos to 

72 \7 / 2 


dg = 4 sin i cos i 


I cos an - sin 2 i cos to (B23b) 
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dj = 4 sin^i cos i cos 


= 2 sin i 


cos 1 


sin 2 i sin - 2^ - sin 2 i^ sin £&> 


d^ = -4 sin sin^i cos i 


(B23c) 

(B23d) 

(B23e) 
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(a) 6C 10 = 10' 5 . 

Figure 2- Sensitivity of range-rate data to variations in gravitational harmonic coefficients for a 46 - 1850 kilometer orbit for which u i 

i 0 = 15°, and M„ = 0°. 
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(d) 5C 22 = 10‘ 7 . 
Figure 2.- Continued. 
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(e) 6C 30 = 10' 5 . 


Figure 2 - Continued. 
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(a) 6C 10 = 10' 5 . 

Figure 3.- Sensitivity of range data to variations in gravitational harmonic coefficients for a 46 - 1850 kilometer orbit for which w, 

in = 15°, and M„ = 0°. 
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(c) 6C 20 = 1(T 6 . 


Figure 3.- Continued. 
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(a) 6C 10 = lCT 5 . 

Figure 6.- Effect of inclination on sensitivity of range-rate data for a 46 - 1850 kilometer orbit for which w 0 = 0°, = 210°, and M 0 = 0°. 


59 






.0 


0 

-.01 

.01 

SA^, 
m/sec 0 


-.01 


.01 

0 


-.01 


64 









"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


